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ABSTRACT 

We report the discovery of two exoplanets transiting high-jitter stars. HAT-P-32b orbits the bright 
V=11.289 late-F-early-G dwarf star GSC 3281-00800, with a period P = 2.150008 ± 0.000001 d. The 
stellar and planetary masses and radii depend on the eccentricity of the system, which is poorly 
constrained due to the high velocity jitter (~ 80 ms -1 ). Assuming a circular orbit, the star has a 
mass of 1.16±0.04M©, and radius of 1.22±O.O2i? , while the planet has a mass of 0.860 ± 0.164 Mj, 
and a radius of 1.789 ± 0.025 When the eccentricity is allowed to vary, the best-fit model has 
e = 0.177 ± 0.079 and results in a planet which is close to filling its Roche Lobe. We determine 
an analytic approximation for the transit-inferred radius of an eccentric planet which fills its Roche 
Lobe; including the constraint that the planet cannot exceed its Roche Lobe results in the following 
best-fit parameters: e = 0.163±0.061, M p = 0.94±0.17Mj, R p = 2.04± O.lOiZj, M* = 1.18±^ 07 M e 
and R± = 1.39 ± 0.07 R Q . The second planet, HAT-P-33b, orbits the bright V=11.188 late-F dwarf 
star GSC 2461-00988, with a period P = 3.474474 ± 0.000001 d. As for HAT-P-32, the stellar and 
planetary masses and radii of HAT-P-33 depend on the eccentricity, which is poorly constrained due 
to the high jitter (~ 50 ms^ 1 ). In this case spectral line bisector spans are significantly anti-correlated 
with the radial velocity residuals, and we are able to use this correlation to reduce the residual rms 
to ~ 35ms _1 . We find the star has a mass of either 1.38 ± 0.04 M or 1.40 ± O.IOMq, and a radius 
of either 1.64 ± 0.03 R Q or 1.78 ± 0.28 R Q , while the planet has a mass of either 0.762 ± 0.101 Mj or 
0.763 ± 0.117 Mj, and a radius of either 1.686 ± 0.045 i?j or 1.827 ± 0.290 for an assumed circular 
orbit or for the best-fit eccentric orbit respectively. Due to the large bisector span variations exhibited 
by both stars we rely on detailed modeling of the photometric light curves to rule out blend scenarios. 
Both planets are among the largest radii transiting planets discovered to date. 

Subject headings: planetary systems — stars: individual ( HAT-P-32, GSC 3281-00800, HAT-P-33, 
GSC 2461-00988 ) techniques: spectroscopic, photometric 



1. INTRODUCTION 

One of the most significant findings that has resulted 
from the study of transiting exoplanets (TEPs) over the 
past decade is the discovery that some close-in "hot 
Jupiters" have radii that are substantially larger than 
what was thought to be theoretically possible. The 
three most inflated TEPs, including WASP-17b (R = 
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1.99 ± 0.08 J?i: IAnderson et^[20Tl . WASP-12b (R = 
1.79 ± 0.09i^; lHe"bb et al.l l200iE~ and TrES-4b (i? = 
1.78 ± 0.09 R j ; ISozzetti et al.ll2009D have radii that are 
as much as 50% larger than exp ected from, for example, 
the coreless lFortnev et al.l ()2007f) models. Recently it has 
become clear that the degree to which TEPs are inflated 
is correlated with the planet equi l ibrium temperature 
dFortnev et al.l 120071: l Enoch et all 12011b Kovacs et alJ 
2010 HFaedi et al.ll201ll:lBekv et al.ll201lULaughlin et al.l 
201 lh and anti-correlated with stel l ar metallicity 
(IGuillot et all 120061: lEortnev et all 120071 : IBurrows et all 
120071: lEnoch' et al. ll201lHBekv et al.l2011f) . Several mech- 
anisms that might explain the correlation with equi- 
librium temperature in p a rticular have been proposed 
dBodenheimer et al.l 120011: IGuillot and Showman] 120021 : 
iBatvgin and Stevenson! I2010h . though to date this issue 
remains unresolved. 

In this work we report the discovery of two new 
TEPs, HAT-P-32b and HAT-P-33b (orbiting the stars 
GSC 3281-00800 and GSC 2461-00988, respectively), 
with radii among the largest found to date. These plan- 
ets have high equilibrium temperatures, supporting the 
aforementioned correlation. Both of these highly inflated 
hot Jupiters were discovered by the H ungarian- made Au- 
tomated Telescope Network fHATNet: [Bakos et al.ll2004D 
survey for TEPs orbiting bright stars (9 < r < 14.5). 
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HATNet operates six wide-field instruments, including 
four at the Fred Lawrence Whipple Observatory (FLWO) 
in Arizona, and two on the roof of the hangar servic- 
ing the Smithsonian Astrophysical Observatory's Sub- 
millimeter Array, in Hawaii. 

Although the two planets presented here were among 
the first candidates identified by HATNet, with discov- 
ery observations dating back to 2004, they proved to be 
difficult to confirm due to the significant radial veloc- 
ity (RV) jitter exhibited by the stellar hosts (78.7 ms -1 
and 55.1ms" 1 for HAT-P-32 and HAT-P-33, respec- 
tively), which limits the power of the traditional spec- 
tral line bisector technique used to rule out blend scenar- 
ios. We argue that the jitter is astrophysical in origin, 
and likely related to convective inhomogeneities which 
vary in time, perhaps due to time-varying photospheric 
magnetic fields, and we conduct detailed blend model- 
ing of the observations to confirm the planetary nature 
of these systems. The high jitter values for both stars 
result in poor constraints on the orbital eccentricities of 
the two systems. While we could assume circular orbits 
as has been done for other TEPs, we choose not to do so 
because several eccentr ic short-period TEPs ha ve been 
disc overed (e.g. XO-3 b, IJohns-Krull et~aT1[200l WASP - 
14b. Uoshi efall [20091 : HAT-P-21b. iBakos et al.ll2010bf ). 
so the possibility that either planet is eccentric should 
be included in the parameter uncertainties. Importantly 
the inferred planetary radii depend strongly on the or- 
bital eccentricities. In particular, if HAT-P-32b has an 
eccentric orbit its radius may be larger than any other 
known TEP. 

The structure of the paper is as follows. In Section [2] 
we summarize the detection of the photometric transit 
signal and the subsequent spectroscopic and photomet- 
ric observations of each star to confirm the planets. In 
Section [3] we analyze the data to rule out false positive 
scenarios, and to determine the stellar and planetary pa- 
rameters. Our conclusions are discussed in Section |U 

2. OBSERVATIONS 

2.1. Photometric detection 

Table [T] summarizes the HATNet discovery observa- 
tions of each new planetary system. The HATNet im- 
ages were processed and reduced to trend- filtered light 
curves following the procedure described bv IBakos et al.l 
(|2010ah . The light curves were searched for periodic 
box -shaped signa l s usin g the Box Least-Squares (BLS; 
see iKovacs et all 12002ft method. For HAT-P-33 sup- 
porting observations of the discovery were o btained with 
the Wise-HAT fWHAT; IShporer et alJfeOOQl ) telescope at 
Wise Observatory in Israel; these were analyzed in the 
same manner as the HATNet images. We detected signif- 
icant signals in the light curves of the stars summarized 
below: 

• HAT-P-32 - GSC 3281-00800 (also known as 
2MASS 02041028+4641162; a = 02 h 04 m 10.27s 
6 = +46°41'16.2"; J2000; V=l 1.289 iDroege et al.l 
2006). A signal was detected for this star with an 
apparent depth of ~20.3mmag, and a period of 
P =2. 1500 days (see Figure HJ left). 

• HAT-P-33 - GSC 2461-00988 (also known as 
2MASS 07324421+3350061; a = 07 h 32 m 44.20s, 



8 = +33°50'06.2"; J2000; V=11.188 IDroege et al.1 
l200l . A signal was detected for this star with 
an apparent depth of ~8.5mmag, and a period of 
P =3.4745 days (see Figure [IJ right). 



2.2. Reconnaissance Spectroscopy 

High-resolution, low-S/N "reconnaissance" spectra 
were obtained for both HAT-P-32 and HAT-P-33 us- 
ing the Harvard- Smithsonian C enter for Ast rophysics 
(CfA) Digital Speedometer (DS; iLathaml 11992ft on the 
FLWO 1.5 m telescope. These observations, which are 
summarized in Table [2] were reduc ed and analyzed fol - 
lowing t he procedure describ ed by iTorres et al.l (2002) 
(see also lLatham et all 120091 ) . We find that both stars 
show no velocity variation at the lkms -1 precision of 
the observations, and all spectra are consistent with 
single, moderately-rotating dwarf stars. For HAT-P- 
32 the stellar atmospheric parameters that we find, as- 
suming solar composition, are T c g+ = 6500 ± 100 K, 
logg* = 4.5 + 0.25, wsini = 21.9 ± l.Okms -1 , and 
7RV = -23.21 ± 0.26kms -1 . For HAT-P-33 we find 
T Gff * = 6500 + 100 K, log = 4.0 + 0.25, vsini = 
15.6+ l.Okms" 1 , and 7 RV = 23.03 ± 0.28 kms -1 . 

In addition to the DS observations for HAT-P-33, 
we also obtained several initial reconnaissance observa- 
tions of this target with the SOPHIE spectrograph on 
the Observatoire de Haute-Provence 1.93 m telescope. 
These observations showed ~ 100 ms -1 scatter, with 
only a very faint hint of phasing with the photometric 
ephemeris, hinting on the possibility that the system is a 
blend. Based on these observations we postponed further 
follow-up of the target for several years. 

2.3. High resolution, high S/N spectroscopy 

We obtained high-res olution, high-S/ N spectra for 
both stars using HIRES (jVogt et al.lll994D on the Keck-I 
telescope in Hawaii. For HAT-P- 32 we gathered 28 spec- 
tra with the L> absorption cell (jMarcv fc Butlerl 11992ft 
between 2007 August and 2010 December, together with 
a single l2-free template spectrum. We rejected one low- 
S/N spectrum for which we were unable to obtain a 
high-precision RV measurement, and exclude from the 
analysis two spectra which were obtained during tran- 
sit and thu s may be affected b y the Rossiter-McLaughlin 
effect (e.g. lOueloz et all 12000ft . For HAT-P-33 we gath- 
ered 22 spectra with the I2 cell between 2008 Septem- 
ber and 2010 December, and two template spectra. The 
HI RES spectra were r educed to barycentric RVs follow- 
ing [Butler]eFaI] (|1996l ). The resulting measurements are 
given in Tables [3] and H for HAT-P-32 and HAT-P-33 
respectively. The phased data, along with our best fit 
models for both circular and eccentric orbits, are dis- 
played in Figures [2] and [3] 

For both candidates the RV residuals from the best-fit 
model show significant scatter greatly in excess of what 
is expected based on the measurement uncertainties. For 
HAT-P-32 the residual RMS is 80.3 ms- 1 or 75.0 ms- 1 
for circular and eccentric models respectively, while the 
RMS expected from instrumental variations plus pho- 
ton noise is only 15.7ms _1 . For HAT-P-33 the residual 
RMS is 55.7ms _1 or 54.1ms -1 , again for circular and 
eccentric models respectively, while the expected RMS is 
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Fig. 1.— Left: unbinned light curve of HAT-P-32 from HATNet (top see the text for details), folded with the period P = 2.1500085 days. 
The solid line shows a transit-model fit to the light curve (Section 13. 4{ . We also show a zoomed-in view of the transit (bottom), with the 
dark points showing the light curve binned in phase with a binsize of 0.002. Right: same as left, here we show HAT-P-33. 



TABLE 1 

Summary of photometric observations 



Instrument /Field 


Datc(s) 


Number of Images' 1 


Mode Cadence (s) 


Filter 


HAT-P-32 










HAT-7/G125 


2005 Sep-2006 Feb 


7077 


330 


/ band 


HAT-8/G125 


2005 Sep-2005 Oct 


816 


330 


/ band 


HAT-6/G126 


2004 Dec-2005 Mar 


1365 


330 


/ band 


HAT-9/G126 


2004 Dec-2005 Mar 


1505 


330 


/ band 


KeplerCam 


2007 Sep 24 


602 


33 


z band 


KeplerCam 


2007 Oct 22 


489 


53 


z band 


KcplcrCam b 


2007 Oct 23 


161 


38 


z band 


KeplerCam 


2007 Nov 06 


759 


28 


2 band 


KeplerCam b 


2007 Nov 18 


551 


38 


2 band 


KeplerCam 


2007 Nov 19 


665 


38 


2 band 


KcplcrCam b 


2007 Dec 03 


817 


33 


2 band 


KeplerCam 


2007 Dec 04 


596 


33 


g band 



HAT-P-33 

HAT-6/G176 

HAT-9/G176 

WHAT/G221 

KeplerCam 

KeplerCam 

KeplerCam 

KeplerCam 

KeplerCam 

KeplerCam 

KeplerCam 



2004 Nov-2005 Oct 
2004 Nov-2005 Oct 
2004 Jan-2005 May 

2006 Dec 29 

2007 Nov 24 

2008 Mar 25 
2008 Nov 23 
2008 Nov 30 
2011 Feb 14 
2011 Feb 21 



2754 
4383 
5439 
413 
274 
260 
668 
257 
346 
336 



330 
330 
330 
33 
57 
29 
28 
32 
57 
58 



band 
band 
band 
band 
band 
band 
band 
band 
band 
band 



a For HATNet and WHAT observations this includes images which were rejected by the photometric reduction 
pipeline. For KeplerCam observations this excludes images which were rejected by the photometric reduction 
pipeline, but includes images which were rejected via cr-clipping during the fitting procedure. 

Out of transit observation used to constrain the presence of a secondary transit event during blend modeling. 



7.9ms' 1 . Because of this high "jitter" we gathered sub- 
stantially more high-S/N spectra for these targets than 
we do for typical HATNet candidates. The reason for 
this is two-fold: the high jitter could have been due to 
additional Keplerian motion due to the presence of ad- 
ditional planets in the systems which would be revealed 
by further observations, and in the presence of high jit- 
ter more observations are needed to precisely determine 
the orbital parameters of the system. We consider fur- 
ther the origin of the jitter for each object in Section 1331 
concluding that it is most likely not due to additional 
planets, but rather due to stellar activity. 

For each spectrum we also calculated the spectral line 
bisector span (BS) following the method described in §5 



of lBakos et al.l (120071) and the S activity inde x calibr ated 
to the scale of iVaughan. Preston fc Wilson! (119781) fol- 
lowing the procedure of llsaacson fc Fischerl (|2010f) . Note 
that for the BS we follow the convention BS = V\ — v% 
where v\ is the velocity of a point on the bisector near 
the continuum level and v 2 is the velocity of a point on 
the bisector near the line core. The BS and S values are 
plotted in Figures [5] and [3J For both stars the BS and 
S values do not phase with the photometric ephemeris; 
however the scatters in the BS values are comparable 
to the RV semiamplitudes, so we are not able to rule 
out blend scenarios for either object based on the bisec- 
tors. Instead we rely on detailed blend-modeling of the 
observations to rule out these scenarios as described in 
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TABLE 2 

DS RECONNAISSANCE SPECTROSCOPY 
OBSERVATIONS 



JD - 2400000 


RV a 


<x R ,v b 




(kms- 1 ) 


(kins" 1 ) 


HAT-P-32 






ooyoo.UUUo 


-24.58 


1.37 




-24.02 


1.05 


54016.7473 


-23.51 


0.75 


54070.7298 


-23.80 


1.08 


54072.7003 


-24.07 


0.71 


54075.8081 


-22.86 


0.62 


54077.7263 


-23.25 


0.75 


54421.8025 


-21.91 


0.81 


TZ A 1 O 707H 

54422.7370 


-23.67 


0.92 


54423.7596 


-23.03 


0.86 


54424.7471 


-23.80 


0.78 


54425. ozlU 


-22.60 


0.77 


54427.7552 


-22.34 


1.00 


54430.6748 


-20.83 


1.70 


54726.8681 


-23.50 


0.90 


HAT-P-33 






53864.6428 


22.10 


0.62 


53865.6274 


21.91 


0.48 


53866.6480 


22.41 


0.48 


53873.6285 


24.07 


0.63 


54041.9844 


23.21 


0.76 


54043.9879 


24.41 


0.79 


54047.0178 


22.75 


0.64 


54047.9772 


23.54 


0.50 


54048.9252 


22.84 


0.66 



a The measured heliocentric RV of the 
target quoted on the native CfA system. 
Our best guess is that 0.14kms — 1 should 
be added to put the CfA velocities onto an 
absolute system defined by observations of 
minor planets. 

The RV measurement uncertainty. 

Section EH For HAT-P-33 we found that the RV residu- 
als are correlated with the BS (see also Section [3T2]) . and 
were thus able to use the BS to correct the RVs, reducing 
the effective jitter to ~ 35 ms -1 . 

2.4. Photometric follow-up observations 

We conducted high-precision photometric observations 
of HAT-P-32 and HAT-P-33 using the KeplerCam CCD 
camera on the FLWO 1.2 m telescope. The observations 
for each target are summarized in Table [TJ For HAT- 
P-32, in addition to observations taken during transit, 
we also obtained three sets of observations taken out of 
transit at the predicted times of secondary eclipse (as- 
suming a circular orbit). These data are excluded from 
the global analysis described in Section 13.41 but are used 
in Section [3. 31 in ruling out blend-scenarios. 

The reduction of the KeplerCam ima ges was performed 
as described by iBakos et al.1 (|2010aD . We performed 
External Parameter Decorrelatio n (EPD) and used the 
Trend Filtering Algorithm (TFA; iKovacs et al"1l2005f ) to 
remove trends simultaneo usly with the light curve mod- 
eling (for more details, see lBakos et al.|[2010af ) . The final 
time series, together with our best-fit transit light curve 
models, are shown in Figure |4j the individual measure- 
ments are reported in Tableland Table [6] for HAT-P-32 
and HAT-P-33 respectively. 

3. ANALYSIS 
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Fig. 2. — Top panel: Kcck/HIRES RV measurements for 
HAT-P-32 shown as a function of orbital phase, along with our 
best-fit circular (solid line) and eccentric (dashed line) mode ls as 
determined from our global modeling procedure (Section CO); see 
Table \E\ Zero phase corresponds to the time of mid-transit. The 
center-of-mass velocity has been subtracted. Second panel: Veloc- 
ity O — C residuals from the best fit circular orbit. Jitter is not 
included in the error bars. Third panel: Velocity O — C residu- 
als from the best fit eccentric orbit. Fourth panel: Bisector spans 
(BS), with the mean value subtracted. The meas urement from the 
template spectrum is included (see Section l3.3ll . Bottom panel: 
Chromospheric activity index S measured from the Keck spectra. 
Note the different vertical scales of the panels. Observations shown 
twice are represented with open symbols. 



3.1. Properties of the parent stars 

Planetary parameters, such as the mass and radius, 
depend strongly on the stellar mass and radius, which 
in turn are constrained by the observed stellar spec- 
tra as well as the light curves and RV cur ves. We fol- 
lowed a n iterative procedure, described by IBakos et al.1 
(|2010al) . to determine the relevant stellar parame- 
ters. The procedure involves iterating between infer- 
ing stellar atmospheric parameters (including the effec- 
tive temperature T c ff*, surface gravity log 5*, metallicity 
[Fe/H], and projected rotation velocity wsini) from the 
Keck/HIRES template spectr um using the Spectroscopy 
Made E asy package (SME: IValenti fc Piskunovl [19961) 
and the IValenti fc Fischer] (|2005l ) atomic line database, 
and modeling the light curves and RV curves (see 
Section 13^4)) to determine the stellar density p+. At a 
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TABLE 3 

Relative radial velocities, bisector spans, and activity index 
measurements of hat-p-32. 



BJD a 


RV b 


CRV 


C 


BS 


C"BS 




s d 


Phase 


(2,454,000+) 


(m s 


i \ 
J 


(m s" 




(m s 


J 


(m s 


-1 \ 
) 










336.95660 


-12 


.33 


13. 


77 


-4 


.23 


22 


.35 


0. 


.240 





.168 


337.12852 


-51. 


.77 


12. 


83 


-35 


.38 


13 


.61 





.238 





.248 








' 


" 


8. 


.77 


12 


.61 


0. 


233 





.621 




91. 


.89 


13. 


58 


22. 


03 


15. 


.06 


0. 


.232 





.625 


339.12681 


-32. 


.01 


12. 


59 


-5. 


.23 


17. 


.10 


0. 


239 


0. 


.177 


oon n o o o rr 


60. 


.33 


13. 


50 


31. 


99 


9 


.65 


0. 


239 


0. 


.547 


i a a ok /i on 


-56. 


.05 


12. 


94 


-7. 


.56 


12 


.17 


0. 


232 


0. 


.469 


o45.14Uoz 


-7. 


.96 


12. 


14 


56. 


66 


19. 


.42 


0. 


235 





.974 


397.81437 


-75. 


.23 


14. 


22 


-36. 


.38 


9 


.63 


0. 


236 


0. 


.474 


398.09647 


35 


.38 


14. 


70 


-41. 


.76 


8 


.10 


0. 


234 





.605 


427.84636 


55. 


.80 


14. 


.45 


-55. 


,18 


12. 


.13 


0. 


235 


0. 


.442 


428.85791 


312. 


16 


15. 


38 


-57. 


.89 


15. 


.15 


0. 


240 


0. 


.912 


,i or\ nn no 


-3. 


.17 


14. 


43 


-29. 


.84 


8 


.73 


0. 


233 


0. 


.407 


4oo.yU2b( 


-11. 


.30 


14. 


12 


-35. 


,11 


17. 


.49 


0. 


231 





.491 


458.93801 


144. 


16 


16. 


.08 


-9. 


.85 


11. 


.93 


0. 


232 


0. 


.903 


460.86208 


103. 


67 


18. 


28 


-9. 


.84 


11 


.05 


0. 


230 





.798 


548.72493 


88. 


.38 


16. 


21 


-19. 


.13 


14. 


.89 


0. 


240 


0. 


.664 


635.11482 


73. 


.67 


16. 


94 


155. 


85 


24. 


.11 


0. 


236 


0. 


.845 


636.11275 


-187. 


30 


17. 


27 


100. 


10 


16 


.18 


0. 


235 


0. 


.310 


724.02617 


-170. 


37 


15. 


27 


7. 


.97 


13. 


.53 


0. 


237 


0. 


.199 


725.07143 


37. 


.66 


16. 


62 


89. 


82 


14. 


.96 


0. 


241 


0. 


.686 


727.11674 


40. 


.79 


16. 


16 


90. 


24 


15. 


.73 


0. 


241 





.637 


777.92375 


-153. 


16 


15. 


.28 


-28. 


.12 


10 


.06 


0. 


232 


0. 


.268 


810.82619 


150. 


61 


17. 


27 


-20. 


.45 


12. 


.98 


0. 


230 


0. 


.571 


838.93902 


46. 


.30 


19. 


26 


-33. 


.20 


13. 


.87 


0. 


225 





.647 


1192.92153 


-143. 


71 


19. 


01 


-28. 


.23 


13. 


.83 


0. 


236 


0. 


.289 


1250.81416 


-211. 


65 


17. 


65 


-30. 


.15 


10. 


.29 


0. 


238 


0. 


.216 


1376.10348 


-80. 


.92 


17. 


65 


6 


.54 


11 


.60 


0. 


241 





.490 


1544.91306° 


-186. 


72 


18. 


93 


-81. 


.83 


21. 


.62 


0. 


224 


0. 


.006 



Note. — For the iodine-free template exposures we do not measure the 
RV but do measure the BS and S index. Such template exposures can be 
distinguished by the missing RV value. 

a Barycentric Julian dates throughout the paper are calculated from Coordi- 
nated Universal Time (UTC). 

The zero-point of these velocities is arbitrary. An overall offset 7 re i fitted to 
these velocities in Scction l3.4l has not been subtracted. 

c Internal errors excluding the component of astrophysical jitter considered in 
Section l3~4l 

^ Chromospheric activity index calibrated to the scale of 

IVaughan, Preston fc Wilson! U978D following Tlsaacson fc Fischerl 1 120101 ). 
c Observation during transit which was excluded from the analysis. 



given cycle in the iteration we combine our esti mates of 
T eS * [Fe/H] and p+ with the Yonsei-Yale (YY; lYi et al.l 
[2001] ) series of stellar evolution models to determine the 
stellar mass, radius, and surface gravity among other pa- 
rameters. If the resulting surface gravity differs signifi- 
cantly from the value determined from the spectrum with 
SME, we repeat the analysis fixing the surface gravity in 
SME to the new value. 

For each star, the initial SME analysis, in which the 
surface gravity was allowed to vary, yielded the following 
values and uncertainties: 

• HAT-P-32: T eff * = 6001 ±88K, [Fe/H] = -0.16 ± 
0.08 dex, log 5* = 4.02 ± 0.07 (cgs), and vsini = 
21 ± O.Skms- 1 . 

• HAT-P-33: T ef f* = 6234±114K, [Fe/H] = -0.04± 
0.08dex, logg* = 3.86 ± 0.10 (cgs), and vsini = 
13.9i0.5kms- 1 . 

As described in Section 13.41 the inferred stellar densi- 
ties, and hence radii, for HAT-P-32 and HAT-P-33 de- 
pend strongly on the orbital eccentricities, which due to 



the high RV jitters are poorly constrained. For each sys- 
tem we conducted separate analyses, first assuming a cir- 
cular orbit and then allowing the eccentricity to vary. In 
each case we obtained a new value of logg* which we 
held fixed during a second SME iteration. The final stel- 
lar parameters for each star, assuming both circular and 
eccentric models, arc listed in Table [3 

The inferred location of each star for both the circular 
and eccentric models, in diagrams of a/i?* (which is re- 
lated to p+) versus T e g+, analogous to the classical H-R 
diagram, is shown in Figures [SJS1 In each case the stellar 
properties and their la and 2a confidence ellipsoids are 
displayed against the backdrop of model isochrones for a 
range of ages, and the appropriate stellar metallicity. 

As a check on our stellar parameter determinations 
we compare the measured photometric colors of the two 
stars to the predicted values based on the models. HAT- 
P-32 has J- K = 0.281±0 OOcQ and V-I c = 0.62±0.10 
(from TASS; lDroege et al"1 l2006), whereas the predicted 

11 Taken from 2MASS HSkrutskie et al.1 12006T) and converted 
to the ESQ ph otometric system using the transformations by 
Carpenter (2001). 
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TABLE 4 

Relative radial velocities, bisector spans, and activity index 
measurements of hat-p-33. 



BJD a RV b <t rv c BS ctbs S d Phase 





(ms 


! 


, _ 
(ms 


J 


(ms 


l\ 
J 


(ms 


> 










728.11380 










-60 


23 


2 


93 





171 





822 


728.12071 


107 


70 


6 


84 


-56 


67 


1 


32 





172 





824 


778.03749 


-88 


19 


6 


47 


-59 


99 


4 


67 





174 





190 


779.11168 


39 


67 


7 


79 


-16 


02 


6 


12 





173 





499 


780.05973 


106 


79 


7 


40 


-2 


11 


1 


05 





173 





772 


791.11863 


76 


24 


7 


97 


-1 


56 


4 


00 





174 





955 


805.96712 


-168 


18 


7 


63 


83 


01 


6 


79 





171 





229 


807.00978 


-76 


71 


9 


16 


74 


87 


7 


06 





170 





529 


807.99835 


125 


68 


6 


92 





18 


3 


79 





172 





813 


809.98993 


13 


78 


9 


01 


-69 


27 


5 


86 





170 





387 


810.88079 


125 


60 


8 


34 


-31 


31 


1 


57 





172 





643 


838.95648 


11 


95 


7 


70 


5 


62 


5 


71 





172 





724 


865.01789 


-25 


12 


8 


25 


16 


88 


4 


31 





171 





224 


954.83936 


-67 


35 


7 


49 


7 


29 


1 


92 





175 





076 


955.85129 


-52 


75 


8 


95 


-4 


96 


7 


02 





173 





367 


1192.01643 










13 


37 


2 


83 





172 





339 


1192.02324 


-93 


30 


7 


19 


22 


62 


3 


16 





173 





341 


1193.09374 


-24 


35 


7 


52 


60 


32 


5 


56 





175 





649 


1193.93089 


-8 


61 


8 


57 


38 


61 


5 


33 





170 





890 


1251.91465 


13 


52 


7 


72 


16 


56 


4 


23 





172 





578 


1468.06245 


22 


33 


7 


86 


91 


25 


6 


82 





173 





789 


1470.09079 


18 


65 


8 


16 


-127 


71 


8 


71 





174 





372 


1545.09930 


26 


89 


8 


37 


-5 


60 


5 


38 





170 





961 


1545.93211 


-31 


18 


7 


91 


7 


80 


4 


89 





165 





201 



Note. — For the iodinc-frcc template exposures we do not measure the RV 
but do measure the BS and S index. Such template exposures can be distin- 
guished by the missing RV value. 

a Barycentric Julian dates throughout the paper arc calculated from Coordi- 
nated Universal Time (UTC). 

The zero-point of these velocities is arbitrary. An overall offset 7 re i fitted to 
these velocities in Scction |3.4| has not been subtracted. 

c Internal errors excluding the component of astrophysical jitter considered in 
Section \3A\ 

Chromospheric activity index calibrated to the scale of 
IVaughan. Preston feW ilson (1978) following]! saacson fc Fischer! 12010) . 



TABLE 5 

High-precision differential photometry of HAT-P-32. 



TABLE 6 

High-precision differential photometry of HAT-P-33. 



BJD 

(2,400,000+) 


Mag a 


CMag 


Mag(orig) b 


Filter 


BJD 

(2,400,000+) 


Mag a 


CMag 


Mag(orig) b Filter 


54368.74140 


-0.00204 


0.00143 


10.54500 


z 


54099.68885 


0.00058 


0.00085 


9.51222 i 


54368.74180 


0.00090 


0.00144 


10.54840 


z 


54099.69425 


0.00071 


0.00083 


9.51224 i 


54368.74217 


0.00417 


0.00143 


10.55120 


z 


54099.69465 


-0.00123 


0.00083 


9.51016 i 


54368.74255 


-0.00158 


0.00143 


10.54590 


z 


54099.69502 


-0.00319 


0.00083 


9.50882 i 


54368.74296 


0.00016 


0.00142 


10.54740 


z 


54099.69541 


0.00122 


0.00083 


9.51279 i 


54368.74335 


0.00079 


0.00143 


10.54810 


z 


54099.69579 


0.00204 


0.00083 


9.51399 i 


54368.74373 


0.00331 


0.00144 


10.55080 


z 


54099.69619 


0.00045 


0.00083 


9.51219 i 


54368.74411 


0.00071 


0.00142 


10.54810 


z 


54099.69659 


-0.00038 


0.00083 


9.51120 i 


54368.74449 


-0.00132 


0.00142 


10.54610 


z 


54099.69697 


-0.00305 


0.00083 


9.50858 i 


54368.74489 


-0.00129 


0.00141 


10.54610 


z 


54099.69735 


-0.00342 


0.00083 


9.50826 i 



Note. — This table is available in a machine-readable form in the 
online journal. A portion is shown here for guidance regarding its form 
and content. 

a The out-of-transit level has been subtracted. These magnitudes have 
been subjected to the EPD and TFA procedures, carried out simulta- 
neously with the transit fit. 

k Raw magnitude values without application of the EPD and TFA 
procedures. 



Note. — This table is available in a machine-readable form in the 
online journal. A portion is shown here for guidance regarding its form 
and content. 

a The out-of-transit level has been subtracted. These magnitudes have 
been subjected to the EPD and TFA procedures, carried out simulta- 
neously with the transit fit. 

k Raw magnitude values without application of the EPD and TFA 
procedures. 



values are J-K = 0.32 + 0.02 and V-I c = 0.574±0.023 
for the circular model, and J — K = 0.34 ± 0.02 and 
V-Ic = 0.592 ± 0.024 for the eccentric model. In both 
cases the model V — Ic colors are consistent with the 
measurements, but the model J — K colors are slightly 



(~ 2a) redder than the measurements. HAT-P-33 has 
J - K = 0.280 ± 0.030 and V - I c = 0.577 ± 0.087 
(also from 2MASS and TASS), and has predicted values 
of J — K = 0.27 ± 0.02 and V - I c = 0.511 ± 0.022 
for the circular model, and J ~ K — 0.28 ± 0.02 and 
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0.0 0.2 0.4 0.6 0.8 1.0 

Phase with respect to T c 

Fig. 3. — Top panel: Keck/HIRES RV measurements for 
HAT-P-33 shown as a function of orbital phase, along with our 
best-fit circular (solid line) and eccentric (dashed line) mode ls as 
determined from our global modeling procedure (Section I3.4II ; see 
Table [9] Zero phase corresponds to the time of mid-transit. The 
center-of-mass velocity has been subtracted. Second panel: Same 
as top panel, here we have subtracted a linear correlation with the 
spectral line bisector spans (BS) from the RV measurements. This 
correlation was determined simultaneously with the fit, and sig- 
nificantly reduces the residual RMS. For the displayed points we 
subtract the correlation determined in the circular orbit fit. Third 
panel: Velocity O — C residuals from the best fit circular orbit in- 
cluding the BS correlation. Jitter is not included in the error bars. 
Fourth panel: Velocity O — C residuals from the best fit eccentric 
orbit including the BS correlation. Fifth panel: BS, with the mean 
value subtracted. Th e me asurement from the template spectrum is 
included (see Section 1 3 - 3D - Bottom panel: Chromospheric activity 
index S measured from the Keck spectra. Note the different verti- 
cal scales of the panels. Observations shown twice are represented 
with open symbols. 



V 



0.522 ± 0.022 for the eccentric model. The 



measured photometric colors of HAT-P-33 are within lo- 
ot the predicted values for both models. 

Neither HAT-P-32 nor HAT-P-33 shows significant 
chromospheric emission in their Ca II H and K line c ores. 
Following the procedure of Isaac son fc Fischerl (|2010f ). we 
find median log i?HK <|Noves et al.lll984h values of —4.62 
and -4.88 for HAT-P-32 and HAT-P-33, respectively. 

3.2. Stellar Jitter 

Both HAT-P-32 and HAT-P-33 exhibit notably high 
scatter in their velocity residuals. Stellar jitter values of 
78.7ms- 1 and 55.1ms- 1 for HAT-P-32 and HAT-P-33 
respectively must be added in quadrature to their veloc- 
ity errors to achieve reduced % 2 values of unity for the 
best-fit circular orbits (see Section I3.4|) . 

A possible cause of the jitter is the presence of one 
or more additional planets in either system, though this 
would not explain the high scatter seen in the BS mea- 
surements. F igure Fp shows the Lomb-Scargle frequency 
spect ra (L-S; lLolriblTl97l^ IScarridllQSl IPress fc Rvbickl 
1989) of the RV observations for both stars. In both cases 
there is a peak in the spectrum at the transit frequency; 
in the case of HAT-P-32 this is the highest peak in the 
spectrum, while in the case of HAT-P-33 the highest peak 
is at a high frequency alias of the transit frequency (if 
we restrict the search to frequencies shorter than 0.49 d^ 1 
the transit frequency is the highest peak). We note that 
in neither case would the planet be detectable from the 
RV data alone. For HAT-P-32, the false alarm proba- 
bility of detecting a peak in the L-S periodogram with a 
height greater than or equal to the measured peak height 
is ~ 8% (this is determined by applying L-S to simulated 
Gaussian white noise RV curves with the same time sam- 
pling as the observations). For HAT-P-33 the false alarm 
probability of the transit peak is 79% (the false alarm 
probability for the highest peak in the spectrum is 30%). 
We also show the L-S spectra of the RV residuals from 
the best-fit circular orbit models, in order to see if there 
is evidence for additional short period planets in the sys- 
tems. The highest peaks in these spectra are at periods 
of 18.104d and 0.8507d for HAT-P-32 and HAT-P-33 re- 
spectively (for each system there is some ambiguity in the 
peak that depends on the frequency range searched and 
whether or not a multiharmonic period search is used). If 
we assume that each system has an additional planet at 
the above stated periods, the RMS of the RV residuals 
decreases to ~ 40 ms" 1 and ~ 30 ms" 1 for HAT-P-32 
and HAT-P-33 respecitvely. However, these peaks have 
corresponding false alarm probabilities of 8% and 8.9%, 
so they are not statistically significant detections. We 
conclude that while either system may host additional 
short period planets, there are not at present enough RV 
observations to make a statistically significant detection. 

Another potential source of RV and BS jitter is vary- 
ing sca ttered moonlight con taminating the spectra. Fol- 
lowing iKovacs et al.1 (|2010t ) we investigated the possi- 
bility that the BS and RV may be affected by vary- 
ing sky contamination, and found no evidence that this 
is the case. Alternatively, RV and BS jitter might be 
caused by variable contamination from a nearby star. 
Based on our KeplerCam observations of HAT-P-32 and 
HAT-P-33, we can rule out contaminating neighbors with 
Ai < 5 mag at a separation greater than 3", or neigh- 
bors with Az < 2 mag at a separation greater than 1" 
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Fig. 4.— Unbinned transit light curves for HAT-P-32 (left) and HAT-P- 33 (right), acquired with KeplerCam at the FLWO 1.2 m 
telescope. The light curves have been EPD and TFA processed, as described in Bakos ct al. (2010a). The dates of th e eve nts are indicated. 
Curves after the first are displaced vertically for clarity. Our best fits from the global modeling described in Section 13.41 are shown by the 
solid lines. Residuals from the fits are displayed at the bottom, in the same order as the top curves. The error bars represent the photon 
and background shot noise, plus the readout noise. 





Effective temperature [K] 



6000 5500 

Effective temperature [K] 



Fig. 5. — Model isochrones from Yi ct al. (2001) assuming a circular orbit (left) and the best-fit eccentric orbit (right) for the measured 
metallicity of HAT-P-32, and ages of 0.5 Gyr and 1-14 Gyr in 1 Gyr steps (left to right in each plot). The adopted values of T e ff A and a/ R+ 
are shown together with their 1<t and 2a confidence ellipsoids. The initial values of T e fj* and a] from the first SME and light curve 
analyses are represented with a triangle. 



of either star. We cannot rule out fainter close neigh- 
bors which could be responsible for at least some of the 
jitter. Constraints on potential neighbors based on mod- 
cling the photometric light curves are further considered 
in Section I3~3l 

Figure E] compares the RV residuals to the BSs for 
HAT-P-32 and HAT-P-33. For HAT-P-33 there is a clear 
anti-correlation between these quantities (the Spe arman 
rank-order correlation test fe.g. lPress et al.lll992f ) yields 
a correlation coefficient of r s — —0.77 with a false-alarm 



probability of 0.04% assuming a circular orbit, while for 
the best-fit eccentric orbit r s = —0.72 with a 0.1% false 
alarm probability). For HAT-P-32 there is a hint of 
an anti-correlation, but it is not statistically significant 
(r s = —0.30 with a false-alarm probability of 14% as- 
suming a circular orbit, and r s — —0.39 with 4.9% false 
alarm probability for the best-fit eccentric orbit). The 
anti-correlation for HAT-P-33 is an indication that the 
measured jitter for this star is dominated by intrinsic stel- 
lar variations, rather than being due to another planet in 
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Fig. 6. — Same as Figure[5] here we show the results for HAT-P-33. In this case the 0.5 Gyr isochrone is not included. 



TABLE 7 

Stellar parameters for HAT-P-32 and HAT-P-33 





HAT-P- 


32 


HAT 


-P-33 




Parameter 


Value a 


Value 


Value a 


Value 


Source 




Circular 


Eccentric 


Circular 


Eccentric 





Spectroscopic properties 

T cff * (K) 

[Fe/H] 

iisini (kms 1 ) . . . 

%ac (kms -1 ) 

ti mic (kms" 1 ) .... 
7rv (kms -1 ) 

Photometric properties 

V (mag) 

V-I c (mag) 

J (mag) 

H (mag) 

K s (mag) 

Derived properties 

A/* (M ) 

R*(R@) 

log g* (cgs) 



My (mag) 

M K (mag.ESO) . . 

Age (Gyr) 

Distance (pc) .... 



6207 ± 88 


6001 ± 88 


6446 ± 88 


6401 ± 88 


SME b 


-0.04 ±0.08 


-0.16 ±0.08 


0.07 ±0.08 


0.05 ± 0.08 


SME 


20.7 ±0.5 


21 ±0.5 


13.7 ±0.5 


13.8 ± 0.5 


SME 


4.69 


4.3 


5.06 


4.99 


SME 


0.85 


0.0 


0.85 


0.85 


SME 


-23.21 ± 0.26 




23.03 ±0.28 




DS C 



11.289 
0.62 ±0.10 
10.251 ±0.022 
10.024 ± 0.022 
9.990 ±0.022 



1.160 ±0.041 
1.219 ±0.016 
4.33 ± 0.01 
1.97 ±0.15 
4.04 ±0.10 
2.77 ±0.04 
2.7 ±0.8 
283 ±5 



1.176" 1 



H0.043 
-0.070 

1.387 ±0.067 
4.22 ±0.04 
2.43 ±0.30 
3.82 ±0.14 
2.50 ±0.11 

d - 8 -o.s 
320 ± 16 



11.188 
0.577 ± 0.087 
10.263 ± 0.021 
10.061 ± 0.024 
10.004 ±0.018 



1.375 ± 0.040 
1.637 ± 0.034 
4.15 ± 0.01 
4.15 ± 0.33 
3.19 ± 0.10 
2.11 ±0.05 
2.3 ±0.3 
387 ±9 



1.403 ±0.096 
1.777 ± 0.280 
4.09 ±0.11 
4.73±1-|J 
3.06 ±0.35 
1.93 ± 0.34 
2.4 ± 0.4 
419 ± 66 



TASS 

TASS 

2MASS 

2MASS 

2MASS 



YY±a/fi*±SME d 

YY±a/R*±SME 

YY+o/R*+SME 

YY+a/Ri+SME 

YY±a/i?*±SME 

YY±a/i?*±SME 

YY+a/Ri+SME 

YY+a/^+SME 



a The eccentricities of both HAT-P-32 and HAT-P-33 are poorly constrained— both planets arc consistent with circular 
orbits but may also have significant eccentricities. Wc list separately the parameters obtained when a circular orbit is 
fixed, and when the eccentricity is allowed to vary. 

k SME — "Spectroscopy Made Easy" package for the analysis of high-resolution spectra JValcnti &c Pis kunovM1996l) ■ 
These parameters rely primarily on SME, but have a small dependence also on the iterative analysis incorporating the 
isochrone search and global modeling of the data, as described in the text. 

c The mean heliocentric velocity as measured by the CfA Digital Speedometer. The uncertainty is the rms of the 
measured velocities of the star rather than the uncertainty on the mean. 

YY+a/i^^+SME — Based on the YY isochroncs ( Yi ct al. 2001), a/R„ as a luminosity indicator, and the SME results. 



the system, or to the instrumental errors being underes- 
timated. We make use of this correlation in Section l3~4l 
to reduce the scatter in the RV residuals for HAT-P-33 
a nd improve the pr ecision of the fitted parameters. 

iSaar et ail (|1998l) argue that the RV jitter of F stars 
is due primarily to convective inhomogeneities (due for 
example to regions of strong magnetic fields which lo- 
cally suppress convection) which vary in time, rather 
than temperature inhomogeneities (due for example to 
cool starspots) which vary in time. The latter may be 



more important for G and K stars. Both sources of jit- 
ter would result in significant BS variations that may 
be correlated with the RV variations, but the lack of 
detectable out-of-transit photometric variations in the 
HATNet light curves of HAT-P-32 and HAT-P-33, and 
the fact the neither star shows particularly high chro- 
mospheric emission in their Ca II H and K line cores, is 
consistent with the expectation that temporally chang- 
ing convective inhomogeneities, rather than temporally 
changing temperature inhomogeneities, are responsible 
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Fig. 7. — Lomb-Scargle frequency spectra for the RVs (top) and residual RVs from the best-fit circular orbit models (bottom) for 
HAT-P-32 (left) and HAT-P-33 (right). In the top panels we mark the transit frequencies, while in the bottom panels we mark the highest 
peaks in the spectra. In all cases the false alarm probability of finding a peak as high the marked peak in a Gaussian white-noise RV curve 
is greater than 8%. 

for the jitter of these two stars. 

Jitter values as high as those found for HAT-P-32 and 
HAT-P-33 are typical of stars with similar spectral types 
and rotation velocities. From t he v sin z-jitter correlation 
measured bv lSaar et all ([2003P ). the expected jitter for an 
F dwarf with vshii = 20kms _1 is ~ 50ms , while 
for an F dwarf with vsmi = 14kms _1 th e expected 
j itter is ~ 30 ms" 1 . The sample used bv iSaar et al.l 
(|2003f ) to determine this correlation includes only a hand- 
ful of stars with wsini > lOkms" 1 , and the scatter 
about the relation is fairly significant — one star with 
vsini ~ 15kms _1 was found to have a jitter in ex- 
cess of 100 ms -1 . The planet hosting star HAT-P-2, 
which has a similar temperature and rotation velocity to 
HAT-P-32 and H AT-P-33 (T e ff* = 6290 ± 60 K, v sini = 
20.8 ± 0.3kms _1 ; IPal et alJl20Tol ) has been reported to 
have a hig h jitter of ~ 60 ms " 1 based on data from Keck 
and Lick (IBakos et al.l [20071 ). A subsequent analysis by 
I Winn et al.l ()2007[ ). using only Keck/HIRES data, found 
a somewhat lower jitter of ~ 30ms _1 , though this is 
still quite a bit higher than most planet-hosting stars 
discovered to date (e.g. the median jitter of the pre- 
viously published HATNet planets is ~ 7ms" 1 ). The 
primary difference between HAT-P-2 and the two sys- 
tems presented here is that the planet HAT-P-2b is sig- 
nificantly more massive than either HAT-P-32 b or HAT- 
P-33b (HAT-P-2b has M p = 9.09 ± 0.24 M j; IPal et al.l 
120101 ): as a result, the RV semiamplitude of HAT-P-2b 
(K = 984 ± 17ms" 1 ) greatly exceeds the jitter, making 
this planet more straightforward to confirm than either 



HAT-P-32b or HAT-P-33b. 

3.3. Excluding blend scenarios 

Both HAT-P-32 and HAT-P-33 exhibit significant 
spectral line bisector span (BS) variations (Figures [2][3j 
the RMS of the BS is 53 ms" 1 and 51 m s" 1 for HAT-P-32 
and HAT-P-33 respectively) . In neither case are the vari- 
ations in phase with the transit ephemeris, as one might 
expect if the observed RV variation were due to a blend 
between an eclipsing binary system and another star. As 
discussed in Section [3~2l these variations are likely due 
to temporally changing convective inhomogeneities, per- 
haps created by variable photospheric magnetic faculae 
as in the Sun, which we suspect are responsible for the 
significant RV jitter seen in both of these stars. Nonethe- 
less the large BS variations, which are comparable to the 
semiamplitudes of the RV signals, prevent us from using 
the BSs to rule out the possibility that either of these 
systems is a blend. 

To rule out blend scenarios we made u se of the 
blendanal program ([Hartman et al.l 1201 lbt see also, 
lHartman et al.ll2011al ) which models the observed light 
curves, stellar atmospheric parameters, and calibrated 
photometric magnitudes using various blended eclips- 
ing binary scenarios as well as scenarios involving 
a transiting planet system potentially blended with 
light from another star. T he program r elies on a 
combination of the Padova ([Girardi et al.l 120001 ) and 
iBaraffe et al.l ([19981 ) stellar evolution models, the Eclips- 
ing B i nary Orbit Prog r am ( Pop per fc Etiel HMO; lEtzel 
[19811: iNeison fc Davisl [1972J, EBOP) as modified by 
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Fig. 8. — RV residuals from the best-fit circular (top) and eccentric (bottom) model orbits vs. BS for HAT-P-32 (left) and HAT-P-33 
(right). There is a hint of an anti-correlation between these quantities f or HA T-P-32, though it is not statistically significant. For HAT-P-33 
the quantities are clearly anti-correlated. See the discussion in Section I3.2I 



ISouthworth et al.l (i'2004a,b) , an d stellar limb darkening 
parameters from IC aretl (120041) 



It is similar to the 



blender program ( Torres et al . 2005) which has bee n 
used to confirm Kepler planets (e.g. iTorres et al.l 1201 1[ ). 
but with a number of technical di fferences which are de- 
scribed bv lHartman et all (|2011bf ). 
For each object we fit 4 classes of models: 

1. A single star with a transiting planet. 

2. A planet transiting one component of a binary star 
system. 

3. A hierarchical triple stellar system. 

4. A blend between a bright stationary star, and a 
fainter, physically unrelated eclipsing binary. 

Initially we assume that the eclipsing components have 
a circular orbit. For both HAT-P-32 and HAT-P-33 we 
find that the class 1 model (a single star with a transiting 
planet), or the class 2 model with a planet- host star that 
is much brighter than its binary star companion, provide 
better fits (lower \ 2 ) than the class 3 and class 4 models. 
To evaluate the statistical significance with which the 
class 3 and class 4 models may be reje cted, we follow 
the Mo nte Carlo procedure described in iHartman et al.l 
(|2011af >. We find that for HAT-P-32 we may reject both 
the class 3 and 4 models with ~ 13c confidence (the best- 
fit class 4 model consists of a group of stars with similar 
parameters to the best- fit class 3 model) , while for HAT- 
P-33 we may reject the class 3 and 4 models with ~ 6c 
confidence (e.g. Figure HJ. 




Fig. 9. — The tr-level at which the class 3 blend model (hier- 
archical triple stellar system) can be rejected for HAT-P-33 as a 
function of the masses of the two largest stars in the system: Mi is 
the mass of the uneclipsed star, and M2 is the mass of the primary 
star in the eclipsing system. Note that the lack of several kms 
RV variations leads to the constraint Mi > M2 . The best-fit model 
is rejected with ~ 6<r confidence. 

For both HAT-P-32 and HAT-P-33 the eclipsing binary 
star blend scenarios (classes 3 and 4) are excluded in part 
due to the lack of an apparent secondary eclipse or out 
of transit variation that are predicted by models capable 
of fitting the observed primary transits. For HAT-P- 
33 the HATNet light curve provides these constraints, 
while for HAT-P-32 three sets of KeplerCam observations 
collected during predicted secondary eclipses (assuming a 
circular orbit) augment the constraints provided by the 
HATNet light curve. Figure [10] shows a few example 
model light curves for the best-fit class 4 model for HAT- 
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P-32 which illustrate this. 

Because the time of secondary eclipse depends on the 
eccentricity, which is poorly constrained by the RV ob- 
servations, we repeat the blend analysis for both systems 
fixing the eccentricities to the best-fit values as deter- 
mined in Section I3~4l While for HAT-P-32 this scenario 
causes the secondary eclipses to not occur during the out 
of transit KeplerCam observations, the eccentric eclips- 
ing binary results in stronger out of transit variations 
which are ruled out by the HATNet data. In this case 
for HAT-P-32 we may reject the class 3 and 4 models 
with > lla confidence. For HAT-P-33 we may reject the 
class 3 and 4 models with > la confidence. 

For the class 2 models we consider two cases, one in 
which the transiting planet orbits the brighter binary 
star component, and the other in which the planet or- 
bits the fainter binary star component. In both cases 
the spectroscopic temperature and the photometric col- 
ors constrain the brighter star to have a mass of > 1 Mq 
for both systems. We find that the case of the planet 
orbiting the fainter component can be rejected outright 
with > 4cr confidence for both HAT-P-32 and HAT-P- 
33. The case of a planet orbiting the brighter binary 
star component amounts to including third light in the 
fit. When the secondary star contributes negligible light 
to the system the model becomes indistinguishable from 
the best-fit case 1 model, such a model cannot be ruled 
out with the available data. Instead we may place an up- 
per limit on the mass of any binary star companion. For 
HAT-P-32 we find that a putative companion must have 
M < 0.5 Mq with 5cr confidence, while for HAT-P-33 a 
companion must have M < 0.55 Mq with 5cr confidence. 
These translate into upper limits on the secondary to 
primary l/-band luminosity ratio of ~ 1% for both HAT- 
P-32 and HAT-P-33. 

Based on the above discussion we conclude that the 
signals detected in the HAT-P-32 and HAT-P-33 light 
curves and RV curves are planetary in nature. 

3.4. Global modeling of the data 

We modeled the HATNet photometry, the follow-up 
photometry, and the high-precision RV mea surements 
using t he procedure describ ed in detail by iPal et all 
(2008); iBakos et al.l (|2010aD . One significant differ- 
ence from our previous pla net discoveries is that we 
use a Mandcl & Agol (2002) transit model to describe 
the HATNet photometry, rather than a simplified no- 
limb-darkening model. This was necessary due to the 
high S/N HATNet detections, especially for HAT-P- 
32. To describe the follow-up light curves we use a 
iMandel fc Agoll (|2002l ) transit model together with the 
simultaneous External-Parameter-Decorellation (EPD) 
and Trend -Filtering Algorithm model of instrumental 
variations (Bakos e t all[20"l0al ). and we use a Keplerian 
orbit to describe the RV curves. For HAT-P-33 we in- 
clude a term aes x BS in the RV model, where aes is 
a free parameter describing the residual RV-BS corre- 
lation, as we found that this significantly reduces the 
RV residuals. For HAT-P-32 we do not include this cor- 
rection because the residual RV-BS correlation is not 
statistically significant. For each planetary system we 
fit two models, one in which the eccentricity is fixed to 
zero, and another in which the eccentricity varies. The 
parameters for each system are listed in Tables [8][9l Nei- 



ther system has a cl e arly n on-zero eccentricity. Using 
the ILucv fc Sweenevl (|197lD test we find that there is 
a non-negligible ~ 3% probability that the eccentricity 
found for HAT-P-32 arises by chance from a circular or- 
bit, while for HAT-P-33 the probability is ~ 20%. 

Because the resulting planets have large radii (partic- 
ularly when the eccentricity is allowed to vary), the as- 
sumption that each planet has a spherical surface which 
lies well within its Roche Lobe may no longer be valid. In 
particular the fact that a planet cannot exceed its Roche 
Lobe places an upper limit on the transit-inferred radius 
of the planet for a given semi-major axis, eccentricity, 
and mass-ratio. Following the procedure described in 
the appendix, we impose the constraint that the planet 
cannot exceed its Roche Lobe in determining the param- 
eters and errors for each system. Figure [TT] shows an 
approximation of the constraints on R p /a and e for each 
system. In both cases the best-fit solutions are below 
this limit, however the uncertainties on the parameters 
(particularly the upper uncertainties on the eccentrici- 
ties and radii) are reduced by imposing the Roche Lobe 
constraint. 

4. DISCUSSION 

We have presented the discovery of two planets, HAT- 
P-32b and HAT-P-33b, which have radii that are among 
the largest measured to date for all transiting exoplan- 
ets. Figure [T^] shows the location of these two planets 
on a mass-radius diagram, while Figure Q2] shows them 
on an equilibrium temperature-radius diagram. For each 
planet the radius determination depends strongly on the 
eccentricity, which is poorly constrained for both systems 
due to the high stellar jitters. If HAT-P-32b has a circu- 
lar orbit, its radius would be ~ 1.8 Rj, but it could be as 
large as ~ 2.0 Rj if the planet is eccentric. HAT-P-33b 
has a slightly smaller radius of ~ 1.7 Rj if it is on a cir- 
cular orbit, but its radius could also be larger (~ 1.8 i?j) 
if it is eccentric. 

There are only four known transiting planets with 
radii that are comparable to HAT-P-32b or HAT-P 
33b. These are W ASP-17b (R = 1.99 ± 0.08 Rj 
Anderson et alll20Tl . WASP-12b (R = 1.79 ± 0.09 Rj 



Hebb et all 12009^ TrES-4b (R = 1.78 ± 0.09 Rj 
Sozzetti et all I2009D. and K epler- 7b (R = 1.614 ± 
0.015 i?T; lDemorv et al.ll201lD . These four planets orbit 
F stars, as do HAT-P-32b and HAT-P-33b. The rela- 
tively high luminosities of these stars result in relatively 
high planet equilibrium temperatures (T oq > 1600 K in 
all cases). The equilibrium temperature in turn is cor- 
related with planet radius (see references in Section [T]). 
Nonetheless, there must be additional factors which 
cause planets like HAT-P-32b, HAT-P-33b and WASP- 
17b to have radii of 1.7-2 R 3 with T cq ~ 1800 K, while 
the m uch hotter planet WASP-18b (2380K; [Hellier et all 
2009) has a much less inflated radius of only 1.165 
One significant difference between WASP-18b and the 
higher radii planets are the planet masses. While HAT- 
P-32b, HAT-P-33b and WASP-17b have sub- Jupiter 
masses, WASP-18b has a mass of ~ 10 Mj. As seen 
in Figure [T31 it appears that there may be a mass depen- 
dence to the r e q-radius relation, with lower mass planets 
being more strongly impacted by high equilibrium tem- 
perature than higher mass planets. A mass-dependence 
of this type has been predicted by several planet infla- 
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TABLE 8 

Orbital and planetary parameters for HAT-P-32b 



Parameter 


Value a 


Value 




Circular 


Eccentric 



2.150009 ± 0.000001 
2454416.14639 ± 0.00009 
0.1292 ±0.0003 
0.0171 ± 0.0002 
5.32 ± 0.22 
17.84 ±0.03 
0.1508 ± 0.0004 



0.012 



+0.012 



0.108 



0.007 
+0.043 
0.044 
7 ±0.6 



Light curve parameters 

P (days) 2.150008 ±0.000001 

T c (BJD) b 2454420.44637 ± 0.00009 

T w (days) b 0.1295 ± 0.0003 

T12 = T 34 (days) b 0.0172 ± 0.0002 

a A'- 6.05l°°4 

C. A'. 17.80 ±6.03 

Rp/R* 0.1508 ±0.0004 

}? 014+ - 014 

6 = acosi/-R* 0.117loo47 

i (deg) 88.9 ± 0.4 

Limb-darkening coefficients c 

ci,i (linear term) 0.2045 

C2,i (quadratic term) 0.3593 

c u z 0.1527 

c 2 ,z 0.3513 

ci,g 0.4460 

c 2 ,g 0.3107 

RV parameters 

K (ms- 1 ) 122.8 ±23.2 

ecos(w) d 0.000 ±0.000 

esin(w) d 0.000 ±0.000 

e 0.000 ±0.000 

lo (deg) 0±0 

RV jitter (ms" 1 ) 78.7 

Secondary eclipse parameters 

T s (BJD) 2454421.521 ± 0.000 

T a 14 0.1295 ±0.0003 

T M 2 0.0172 ±0.0002 

Planetary parameters 

M p (Mj) 0.860 ±0.164 

R p (Rj) 1.789 ±0.025 

C{Mp,R p ) e 0.10 

p p (gem" 3 ) 0.19 ±0.04 

log 9p (cgs) 2.82t°;^ 

a (AU) 0.0343 ± 0.0004 

T eq (K) 1786 ± 26 

f 0.028 ±0.005 

(F) (10 9 ergs" 1 cm- 2 ) s 2.29 ±0.13 



0.2098 
0.3562 
0.1580 
0.3476 
0.4564 
0.3027 



136.1 ± 23.8 
0.099 ± 0.080 
0.124 ±0.037 
0.163 ±0.061 
52 ± 29 
73.3 



2454417.357 ± 0.109 
0.1653 ± 0.0120 
0.0221 ± 0.0017 



0.941 ±0.166 
2.037 ± 0.099 
0.27 
14+ 003 

2.75 ±0.07 

0344+ - 0004 
"•"J«_0 .0007 

1888 ± 51 
0.027 ± 0.004 
2.86 ± 0.31 



a The eccentricity of HAT-P-32 is poorly constrained. We list separately the parameters ob- 
tained when a circular orbit is fixed, and when the eccentricity is allowed to vary. 
b T c : Reference epoch of mid transit that minimizes the correlation with the orbital period. 
T14: total transit duration, time between first to last contact; T12 — T34: ingress/egress time, 
time between first and second, or third and fourth contact. 

c Values for a quadratic law, adopted from the tabulations by Claret ( 2004) according to the 
spectroscopic (SME) parameters listed in Table [71 

d Lagrangian orbital parameters derived from the global modeling, and primarily determined 
by the RV data. 

c Correlation coefficient between the planetary mass M p and radius R p . 

f The Safronov num ber is given by 6 = 5 (V OS c/V orb ) 2 = (a/ R P )(M P /M+) (sec 
IHansen fc' Barman 200TJ). 

^ Incoming flux per unit surface area, averaged over the orbit. 



14 



Hartman et al. 



Background EB Blend 



Transiting Planet No Blend 




0.025 
0.03 




0.06 



ID 

-0.06 -0.04 -0.02 0.02 0.04 0.06 
0.01 




0.035 



0.005 
0.01 
0.015 
0.02 
0.025 
0.03 
0.035 



3 U.VOO 

0.44 0.46 0.48 0.5 0.52 0.54 0.56 0.44 0.46 0.48 0.5 0.52 0.54 0.56 





-0.4 -0.2 0.2 0.4 



-0.4 -0.2 0.2 0.4 



Phase 



Phase 



Fig. 10. — Comparison of the best-fit class 4 (background eclipsing binary; left) and class 1 (single star with a transiting planet; right) 
blend-models for HAT-P-32. The results are shown for three illustrative light curves: a KeplerCam z band primary transit from 2007 
September 24 (top), a KeplerCam z band out-of-transit light curve from 2007 December 3 (middle), and the Iq band HATNet field G125 
light curve (bottom). The HATNet light curve is folded and binned in phase, using a binsize of 0.002 (note that this is only for display 
purposes, we do not bin the data in the modeling). In each panel the light curve is shown at top together with the model, and the residual 
is shown below. Note that for the KeplerCam observations we plot the EPD/TFA-corrected light curve, this correction is determined 
simultaneously with the fit, as a result there arc slight differences in the plotted KeplerCam light curves for the two classes of blend models. 
The class 4 blend model provides a notably poorer fit to the light curves than the class 1 modcl-this includes slight differences in the 
ingress/egress of the primary transit, and a secondary transit and out of transit variation which arc not seen in the KeplerCam or HATNet 
data. 
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TABLE 9 

Orbital and planetary parameters for HAT-P-33b 



Parameter 


Value 


Value 




Circular 


Eccentric 



3.474474 ± 0.000001 
2455100.50255 ± 0.00023 
0.1836 ±0.0007 
0.0194 ± 0.0002 
6.08 4 



3+0.98 
-0.72 

12.17 ±0.05 
0.1057 ±0.0011 



106+ - 001 

U.1UD_q 0Q1 
n oncr+0.002 

u.ozo_ 002 
7+0.8 
-1.2 



86.7^ 



Light curve parameters 

P (days) 3.474474 ± 0.000001 

T c (BJD) b 2455110.92595 ± 0.00022 

T14 (days) b 0.1839 ± 0.0005 

T12 = T 34 (days) b 0.0195 ± 0.0002 

a /,'. 6.56±° ™ 

C/-R* 12.16 ±0.03 

R p /R* 0.1058 ±0.0011 

0.106™ 

b = acosi/K* 0.325t° °o2 

i (deg) 87.21™ 

Limb-darkening coefficients c 

ci,i (linear term) 0.1762 

C2,i (quadratic term) 0.3768 

ci,z 0.1260 

c 2 ,z 0.3671 

ci,g 0.4149 

c 2 ,g 0.3327 

RV parameters 

K (ms' 1 ) 82. 7± 10.8 

ecos(w) d 0.000 ±0.000 

esin(tj) d 0.000 ±0.000 

e 0.000 ±0.000 

to (deg) 0±0 

RV jitter (ms" 1 ) 34.4 

o B S c -0.814 ±0.164 

Secondary eclipse parameters 

T s (BJD) 2455112.663 ± 0.000 

T s 14 0.1839 ±0.0005 

T M2 0.0195 ±0.0002 

Planetary parameters 

M p (Mj) 0.762 ±0.101 

R p (Rj) 1.686 ±0.045 

C(Mp,R p ) f 0.10 

p p (gcm~ 3 ) 0.20 ±0.03 

logg p (cgs) 2.82 ±0.06 

a (AU) 0.0499 ± 0.0005 

T eq (K) 1782 ± 28 

O g 0.033 ±0.004 

(F) (10 9 ergs _1 cm -2 ) h .... 2.27±0.14 



0.1799 
0.3748 
0.1294 
0.3656 
0.4216 
0.3278 



82.8 ± 12.0 
0.040 ± 0.078 
0.073 ± 0.138 
0.148 ±0.081 
96 ± 119 
36.0 
-0.794 ±0.179 



2455102.330 ± 0.175 
0.2090 ± 0.0480 
0.0230 ± 0.0085 



0.763 ±0.117 

1.827 ± 0.290 

0.34 

15+ - 11 
U - 1O -0.05 

2.75 ±0.13 

0.0503 ± 0.0011 

1838 ± 133 

n30+ - 007 
u.uju_ 005 



2.58 



+0.93 



a The eccentricity of HAT-P-33 is poorly constrained. We list separately the parameters ob- 
tained when a circular orbit is fixed, and when the eccentricity is allowed to vary. 
b T c : Reference epoch of mid transit that minimizes the correlation with the orbital period. 
T14: total transit duration, time between first to last contact; T12 — T34: ingress/egress time, 
time between first and second, or third and fourth contact. 

c Values for a quadratic law, adopted from the tabulations by Claret ( 2004) according to the 
spectroscopic (SME) parameters listed in Table ff] 

d Lagrangian orbital parameters derived from the global modeling, and primarily determined 
by the RV data. 

e Parameter describing a linear dependence of the RVs on the BS values. 
^ Correlation coefficient between the planetary mass M p and radius R p . 

g The Safronov num ber is given by O = 3 (Vo B c/V orb ) 2 = {a/ R P )(M P / Ah) (sec 
Hansen tc Banna n 2007). 
h 



Incoming flux per unit surface area, averaged over the orbit. 
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Fig. 11. — Planet radius normalized to the semi-major axis vs. eccentricity for HAT-P-32b (left) and HAT-P-33b (right). Each point 
corresponds to a single MCMC parameter realization. The solid lines show the maximum allowed transit-inferred planet radius as a function 
of eccentricity corresponding to a planet which fills its Roche-Lobe at periastron (Section [AJ. These are calculated using the median mass- 
ratio and argument of periastron for each system. The wiggles seen for HAT-P-32b are due to numerical noise in the calculation. 
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tion mechanisms (e .g. Fig. 10 of[Guillot 200]| see also 
iBatvein et aI|[20TH ). 

Due to the high jitters of the two stars studied in this 
paper, further high-precision RV observations will not 
significantly constrain the eccentricities of the planetary 
systems. A more promising method would be to ob- 
serve the planetary occultations with the Spitzer space 
telescope. Fortunately both stars are relatively bright, 
(Ks ~ 10.0 mag in both cases), and have expected occul- 
tations deeper than 0.1% in both the 3.5 /mi and 4.6 /iin 
bandpasses. Thus we expect that it should be possible 
to obtain high S/N occultation events (S/N> 10) with 
Spitzer for both systems. 

From standard tidal theory we expect the circulariza- 
tion time-scales to be much shorter tha n the > 2 Gyr ages 
of the systems. Using equation 1 of iMatsumura et al.1 
(2008), and assuming planetary and stellar tidal damp- 
ing factors of Q' p = = 10 6 , the expected circular- 
ization timescales are ~ 3Myr and ~ 30Myr for HAT- 
P-32 b and HAT-P-33b respe ctively. We note, however, 
that iPenev &: Sasselovl (|2011h have recently argued that 
standard tidal theory, which is calibrated from observa- 
tions of binary stars, significantly overestimates the tidal 
interaction between planets and stars and thereby under- 
estimates the true circularization time-scale. It is possi- 
ble that a short period planet may maintain an eccentric 
orbit for the entire life of the system. 

Finally we note that the difficulty of confirming the 
planets presented here illustrates the selection bias im- 
posed on transit surveys by the need for RV confirmation. 
For both planets the HATNet transit detection was clear 
and robust. This, together with the high S/N box-shaped 
transits observed with KeplerCam, motivated us to con- 



tinue obtaining high-precision RVs for these objects de- 
spite the initial RVs not phasing with the photometric 
ephemerides, and the stars showing significant spectral- 
line bisector span variations. If either target had a shal- 
lower, less obviously planet-like transit, it is likely that 
we would not have continued the intensive RV monitor- 
ing necessary for confirmation, and it is also likely that 
we would not have been able to conclusively rule out 
blends based on analyzing the light curves. Moreover, if 
either planet had a significantly lower mass (Saturn-mass 
or smaller), such that the orbital variation could not be 
detected, it is also likely that the planet would not have 
been confirmed. 
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APPENDIX 



CALCULATING THE TRANSIT-INFERRED RADIUS OF AN ECCENTRIC, ROCHE-LOBE FILLING PLANET 

The condition that the surface of a planet cannot extend beyond its Roche-Lobe (assuming the system is not an 
overcontact binary, which is true for the systems considered in this paper since in both cases the inferred stellar radius 
is well within the stellar Roche Lobe for eccentricities that are consistent with the RV curves) sets a maximum limit 
on its size for a given semi-major axis, eccentricity, and star-planet mass ratio. This in turn places an upper limit on 
the radius o f a plan e t infe rred from a transit measurement. Here we briefly review how to calculate this constraint. 

Following IWilsonl (|1979f ) , the binary potential for the general case of nonsynchronously rotating components on a 
non-circular orbit is given by: 



(D 2 



2r\D) 



-1/2 



rX/D 1 



(Al) 



The polar coordinates r, 0, and <p (0 is t ne polar angle) have an origin at the center of one of the binary components 
(in our case we choose the planet), distance is measured in units of the semi- major axis of the relative orbit, A and 
v are direction cosines (A = sin 6 cos <f>, v — cos 6*), D is the instantaneous separation between the planet and star 
(D = 1 — ecosE where e is the eccentricity and E is the eccentric anomaly), q is the mass ratio (q = M+/M p in our 
case), and F is the synchronicity parameter equal to the ratio of the angular rotation velocity of the component at 
the origin (the planet in our case) to the "average" angular velocity of the orbit (2tt/P where P is the orbital period). 
For an eccentric system the tidal interaction drives the components toward s pseudo-s ynchronous rotation, which is 
between the average angular velocity and the angular velocity at periastron (|Hutlll98lT eq. 42): 



F = 



l + fe 2 + fe 4 



-L P 6 

16 c 



(l + 3e 2 + fe 4 )(l-e 2 ) 3 / 2 ' 



(A2) 



We assume pseudo-synchronous rotation for the planet. 

The surface of a gas giant planet is expected to follow a surface of constant potentia l . For the eccentric case the 
value and shape of the surface potential varies with the orbital phase, in this case IWilsonl (fl97l) argues that to a good 
approximation the volume of the object is constant over the orbit, and suggests a procedure, which we adopt, for finding 
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Fig. 12. — Mass-radius diagram of TEPs. HAT-P-32b and HAT-P-33b are indicated. The triangles indicate the parameters for assumed 
circular orbits, while the squares indicate the parameters when the eccentricity is allowed to vary. Filled circles are all other TEPS, and 
filled squares are solar system planets. We also show lines of con stant density (dotted lines running from lower left to upper right) and 
theoretical planet mass-radius relations from Fortncv ct al. (20Q3). 

the surface of a Roche-Lobe filling object at different orbital phases. Because the volume of the Roche Lobe is smallest 
at periastron, it follows that an eccentric planet which fills its Roche Lobe will do so only at periastron. The potential 
fio corresponding to the Roche Lobe of t he planet at periastron can be determined by taking D = 1 — e and finding r 
between and D that minimizes Eq. lAll for v = 0, A = 1. The volume Vo of this surface can be calculated numerically. 
We perform a Monte Carlo integration randomly generating points uniformly distributed within a spherical shell with 
an inner radius ro that satisfies the condition il(ro,9, A) > fin for all 9 and A and an outer radius r\ that satisfies 
fi(ri), 9, A) < fio for all 9 and A, and taking the volume to be equal to 

Vo = \*{rl + /(r? - rg)) (A3) 

where / is the fraction of generated points with O > f^o- The surface potential during transit tt\ may then be 
determined by setting D equal to the appropriate value at transit and finding il such that V(ft = fii) = Vq. We solve 
this using a simple bisection search noting that the volume of the Roche Lobe potential at transit phase (which can 
be determined as at periastron) is greater than Vq. Finally the radius of the planet which would be inferred from a 
transit observation (corresponding to the radius of a circle with area equal to the projected area of the planet surface 
potential as viewed by an Earth-bound observer) is given approximately by 

R p ,rl = y/yz (A4) 

where y = r such that ft(r, 9 = ir/2,<fr = n/2) = fix, z = r such that Vt(r,9 = 0, <p = 0) = fix, and we assume 
for simplicity an edge-on orbit (which is reasonable for a transiting planet). We have also conducted a Monte Carlo 
integration to calculate the projected area of the planet surface potential accounting for inclination, and found the 
difference from the approximation to be negligible (less than 0.1%). 

In principle one may apply the Roche-Lobe constraint in the global modeling of the photometry and RVs by cal- 
culating Rp.RL for each set of trial MCMC parameters, and rejecting the trial if it yields R p /a > R P: rl- In practice 
this is unwieldy because of the slow numerical integrations required by the above procedure. We therefore searched 
for an analytic approximation to R p .rl that depends on the mass ratio q, the eccentricity e and the argument of 
periastron lu. We numerically determined R P: rl over a grid of parameters spanning 10 < q < 10 5 , < e < 0.5, and 
7r/2 < w < 37r/2, and find that the following empirically chosen analytic function reproduces R p ,rl to ~ 1% accuracy 
over this parameter range: 



Rp,RL(q,u,e) w f 1 (q)g 1 (uj)e 2 + / 2 (g)3 2 (w)e + h{q)gz{u) 



(A5) 
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Fig. 13. — Top: T cq -radius diagram of TEPs with 0.1 Mj < M p . HAT-P-32b and HAT-P-33b are indicated with triangles (circular 
models) and squares (eccentric models). The shading of the planets indicates their masses. HAT-P-32b and HAT-P-33b continue the trend 
of higher T cq planets having larger radii. Bottom: the same as above, here we only show planets discovered by the HATNet survey. The 
T eq — radius correlation is clearly apparent when restricted to this sample of planets, which we take as evidence that the correlation is not 
due to combining planets found by different surveys that have different selection biases. 
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TABLE 10 
Best fit parameters for equation IA5I 



PcH'cimctcr 


Value 


(41 


0.0033 ±0.0011 


2 


-0.385 ±0.016 


«3 


-1.673 ± 0.056 


a 4 


-0.487 ±0.036 


as 


-0.00222 ± 0.00032 


06 


-0.2927 ±0.0046 


07 


-1.199 ±0.017 


as 


0.1397 ±0.0086 


ag 


-1.434 x 10" 3 ± 6.9 x 10~ 5 


O10 


-0.30528 ± 0.00097 


on 


-0.9238 ±0.0031 


a 12 


-0.00399 ± 0.00088 



where 



Hi) 


= - exp(ai ln(g) 2 + a 2 ln(q) 


+ a 3 ) 




= 04 cos(a; — 7r/2) + 1 




f2(q) 


= - cxp(a 5 ln(a) 2 + a 6 ln(g) 


+ a 7 ) 


.92 M 


= ag cos(w — 7r/2) + 1 




h(q) 


= exp(a 9 ln(g) 2 + a 10 ln(g) + 


an) 


53 (w) 


= ai2 cos(cj — 7r/2) + 1 





and ai through a 12 are fitted parameters given in Table [TU1 
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